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In this paper, we study the interaction between the two-level atom and a bimodal 
cavity field, namely, two-mode Jaynes-Cummings model when the atom and the 
modes are initially in the atomic superposition state and two-mode squeezed vacuum 
state, respectively. For this system we investigate the atomic inversion, linear entropy 



and atomic Wehrl entropy. We show that there is a connection between all these 



quantities. Also we prove that the atomic Wehrl entropy exhibits behaviors similar 



to those of the linear entropy and the von Neumann entropy. Moreover, we show 



that the bipartite exhibits periodical disentanglement and derive the explicit forms 
of the states of the atom and the modes at these values of the interaction times. 



o 

PACS numbers: 42-50. Dv Jaynes-Cummings model, atomic inversion, two-mode squeezed 
vacuum states, Wehrl entropy, linear entropy. 
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I. INTRODUCTION 

X 

One of the few exactly solvable models in quantum optics, which describes the interaction 
between the radiation field and the matter, is the JaynesCummings model (JCM) [lj. The 
JCM has become experimentally realizable with the Rydberg atoms in the high-Q microwave 
cavities (e.g., see [2j). Various extensions for the JCM have been performed and investigated 
in greater details such as multiphoton 0,0], intensity-dependent (Bf, multimode, e.g., p, 0, 



|8|, ]9|, [lOj, multilevel atoms [ll[] and multiatom interactions [121 ] . One of these extensions 
is the two-mode JCMfTJCM), which has taken a considerable interest in the literatures, 
e.g. 0, 0, @, [Hi EH, Ell- The revival-collapse phenomenon (RCP) occurring in the atomic 
inversion of the TJCM is rather complicated compared to that of the standard JCM (, i.e. a 
single mode interacting with a single atom) in the sense that the revival series is compact and 
each revival is followed with a secondary revival. Furthermore, the RCP associated with the 
TJCM is independent of the initial intensities of the modes. Such behavior has been partially 
explained in [9j, however, an investigation for the occurrence of the secondary revivals is 
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given m [tj||. The quantum phase properties for the TJCM with the initial Schrodinger-cat 
states have been investigated in [3] showing that the phase variances of the single-mode case 
can exhibit RCP about the long-time behavior but with the interaction time several times 
smaller than that of the standard JCM. Furthermore, there is a clear relationship between 
the RCP occurring in the atomic inversion and the behavior of the phase distribution of 
both the single-mode and the two-mode cases. Moreover, for the TJCM it has been shown 
that there is a class of states for which the quadratures squeezing exhibit RCP similar to 
that involved in the corresponding atomic inversion (l3| . 

Quantum information processing provides different way for manipulating information 
than the classical one. This is related to entanglement^which plays an essential role in the 
quantum information 15] such as quantum compu ting (l6| . teleportation [l7j], cryptographic 



181 ] . dense coding [19| and entanglement swapping [20]. Thus intensive efforts have been done 



to understand theoretically and experimentally the entanglement in the quantum systems. 
For instance, the entanglement between two qubits in an arbitrary pure state has been 
quantified by the concurrence [2lJ and Peres-Horodecki measure 22], however, that of the 
mixed states quantified by the infimum of the average concurrence over all possible pure state 
ensemble decomposition. Additionally, the entropic relations are used in investigating the 
entanglement in the quantum system. In this regard von Neuman entropy (NE) [231 ]. linear 
entropy (LE) and the Shannon information entropy (SE]_j24( have been frequently used in 
treating entanglement in the quantum systems. The NE [25| and the LE [3[ have been applied 
to the JCM. It is worth mentioning that the SE involves only the diagonal elements of the 
density matrix and in some cases it can give information similar to that obtained from the 
NE and LE. On the other hand, there is an additional entropic relation, namely, the Wehrl 
entropy (WHE) (26[. This relation has been successfully applied in description of different 
properties of the quantum optical fields such as phase-space uncertainty 0, 28|, quantum 



interference 



decoherence 



, [30| etc. The WHE is more sensitive in distinguishing 

states than the NE since WHE is a state dependent [3l| . The concept of the Wehrl phase 
distribution has been developed and shown that it serves as a measure of both noise (phase- 
space uncertainty) and phase randomization [3lj]. Furthermore, the WHE has been applied 
to the dynamical systems, too. In this respect the time evolution of the WHE for the Kerr- 
like medium has been discussed in 32] showing that the WHE gives a clear signature for the 
formation of finite superpositions of coherent states (cat-like states) as well as the number 
of coherent components taking part in the superposition. For the JCM the WHE gives an 
information on the splitting of the Q-function in the course of the collapse region of the 



atomic inversion as well as on the atomic inversion itself [30, 33[. In the present paper we 
show that the atomic Wehrl entropy (AHE) [33| can be used to quantify entanglement in 
the TJCM when the modes are initially prepared in the maximally entangled states such as 
the two- mode squeezed vacuum state (TMS) [34|. The TMS can be expressed as 



(1) 



n=0 



where C„ 



1 



cosh r ' 



tanh r. We should stress that the TMSs are important states since 



they contain quantum correlations between the different modes that make up the field. Based 
on this fact the TMSs have been used in the continuous- variable teleportation jlBj, quantum 
key distribution j36|, verification of EPR correlations 37 1. etc. Most importantly the TMSs 
can be experimentally generated via the optical parametric oscillator [38(, where squeezing 
can be detected with high efficiency homodyne detector 39]. Now for the TJCM with the 
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initial TMS we investigate the atomic inversion, LE, Wehrl atomic density and AHE. This 
work is motivated by the importance of the TJCM and TMS in quantum optics as well as 
in the quantum information, as we mentioned above. We obtain many of interesting results 
such as the atomic Wehrl density (entropy) can give qualitative (quantitative) information 
on the entanglement in the TJCM. In other words, the AWE can be used as a new measure 
for quantifying entanglement in the quantum systems. Furthermore, we show that the 
bipartite (, i.e. two-mode and atom) can periodically exhibit instantaneous (long-lived) 
disentanglement. Also we derive the explicit forms for the states of the atom and the two- 
mode at these values of the interaction times. The paper is prepared in the following order: 
In section 2 we derive the wavefunction of the system and evaluate the expectation values 
of the atomic variables. In section 3 we investigate the atomic inversion and the LE. In 
section 4 we discuss the atomic Wehrl density and the atomic Wehrl entropy. In section 5 
we summarize the main results. 



II. WAVEFUNCTION OF THE SYSTEM 

In this section we derive the dynamical wavefunction for the system, which consists of 
two-mode field interacting with the two-level atom. Also, we evaluate the expectation values 
for the atomic operators which will be frequently used in the paper. 

The Hamiltonian describing the interaction between two-mode field and the two-level 
atom in the rotating wave approximation takes the form [4(| : 



where 



H = h(H + H l ), (2) 



2 

(3) 



Hi = xia^df + a^al) + A dia 2 S + + a\a\S. 

where ctj (aj) and uij are the annihilation (creation) operator and frequency designated the 

jth mode, the set {S z , S + , S-} is the usual Pauli spin operators, A is the coupling constant 
between the atom and the modes, x is the dispersive part of the third-order nonlinearity of 
the Kerr-like medium, which is assumed to be the same for the two modes. Also uo a denotes 
the atomic transition frequency. 

Throughout the investigation we consider that u a = u>i + u>2 (i.e. the resonance case), 
the modes are initially prepared in the TMS (pQ) and the atom is in the atomic superposition 
state having the form: 

|vl/ a (0)) = cos(^)|T) + sin(^)exp(^)|I), (4) 

where ||) and ||) deonte excited and ground atomic states, respectively; 9 and are phases. 
It is worth reminding that preparing the atom in the coherent superposition states is im- 
portant as a result of its applications in noise quenching by correlated spontaneous emission 



4l| . quantum beats [42j and noise- free amplification [431] • Now under these conditions the 
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dynamical wavefunction for the Hamiltonian system ([3]) can be evaluated in the standard 
way as 

oo 

I ¥(T)) = ^exp(-2^Tn 2 ) [Fi (n, T) || ,n,n) + F 2 (n,T) |j, n + 1, n + 1)] , (5) 

n=0 

where 

Fi(n,T) = C n cos f cos(Tfi n ) + §^ [2r/racos| - (n+l)exp(i</>) sin §] sin(Tfi n ), 

(6) 

F 2 (n,T) = C„exp(i0) sin | cos(7Y2„) — ^ [2r/nexp(z0) sin | + (n + l)cos|] sin(Tf2„) 
and 

fi n = VVn 2 + (n + l) 2 , 7}=$, T = tX. (7) 

A 

All information about the system is involved in the wavefunction (J5J) or in the total density 
matrix p(T) =\ ^/(T))(^/(T)\. As we mainly deal with the atomic subsystem we evaluate 
the atomic reduced density matrix p a (T) via the relation: 

p a (T) = Tr fp(T), (8) 

where the subscript / means that the trace is taken over the field. We close this section by 
evaluating the expectation values for the atomic operators {S z , S x , S y } for (JSJ) as 

oo 

(S X (T)} = \ Px (T) = Re £ exp[2zyr(2n + l)]F?(n + l,T)F 2 (n,T), 

71=0 
OO 

(S y (T)) = \p y {T) =Im E exp[2zryT(2n+ l)]F*(n + l,T)F 2 (n,T), (9) 

n=0 

oo 

(S Z (T)) = \ Pz (T) = \ £ [li^T)! 2 - |F 2 (n,T)| 2 ], 

n=0 

where ^ = 1(5*- + S + ),S y = ^(S + — S_). We use the results obtained in this section 
to make a comparative study between the atomic inversion, LE and AWE in the following 
sections. 



III. ATOMIC INVERSION AND LINEAR ENTROPY 

In this section we investigate the atomic inversion and the LE for the system under 
consideration. We start the investigation with the atomic inversion, which is plotted in Fig. 
1 versus the scaled time T for the given values of the interaction parameters. From this figure 
the obvious remark is that the (S Z (T)) is periodic with period it, regardless of the type of 
the initial atomic state. The origin in this is related to the nature of the TMS in which the 
information is equally distributed among the two modes Furthermore, such periodicity 
is quite similar to that of the two-photon JCM [4] and intensity-dependent JCM [5|. For 
the superposition atomic state (9, (f>) = (ir/2,ir/2), i.e. the dashed curve, one can observe 
that (S Z (T)) oscillates around zero with extreme values ±0.15. Nevertheless, we have found 
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FIG. 1: The (S Z (T)) versus scaled time T for sinh 2 r = 10 and (9,(j),r]) = (0,0,0) (solid curve) 
and (-7r/2, 7r/6, 0) (dashed curve). 



that the system can exhibit atomic trapping, i.e. (S Z (T)} = 0, for (#,</>) = (7r/2, 0), as we 
show shortly. All these facts can be analytically realized for rj = 0. In this case (S Z (T)) in 
(P) can be easily modified as: 



(S Z (T)) 



[cos(2T) - z 2 } cos9 + z sin 6 sin(0) sin(2T) 
2[(1 - z 2 cos(2T)) 2 + z A sin 2 (2T)] cosh 2 r ' 



(10) 



From ([TO]) the smallest value of the atomic inversion for 9 = is (S Z (T)) = — 1/(2 + 2n), 
where n = 2 sinh 2 r, which occurs at T = n/2 (, i.e. collapse time). This means that the 
greater the initial mean-photon number of the field the closer the value of the (S Z (T)) to 
zero at the collapse time. Also it is obvious that for (9,<p) = (tt/2,0) the expression ffTUj) 
tends to zero. We have found that the behavior of the (S Z (T)) is insensitive of the value of 
t]. Finally the information obtained here is important in investigating the entanglement in 
the bipartite, as we show below. 

Now we quantify the entanglement between the two-mode system and the atomic system 
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FIG. 2: The f(T) versus scaled time T for sinh 2 r = 10 and (p,<f>,rj) = (0,0,0) (a), (vr/2, vr/6, 0) 
(b) and (0,0,0.03) (c). 



via the linear entropy, which is defined by 

£(T) = l-Trp*(T). 



(11) 



As a result of the symmetry in the system the linear entropies for the atomic and field 
sub-systems are equivalent. Furthermore, as the system is in an over all pure state, the LE 
is a relevant measure for quantifying entanglement in the bipartite (45|. The £(T) ranges 
from for pure states and/or disentangled states to 0.5 for a maximally entangled states 
For the future purpose, the relation (TTTj) can be expressed in terms of the expectation 
values of the atomic variables (191) as 



pl(T) - P 2 y (T) - pl(T)] . 



(12) 



In Figs. 2 we plot £ (T) versus the scaled time T for given values of the interaction parameters. 
From Fig. 2(a) and Fig. 2(b), where i] = 0, the bipartite is periodically disentangled with 
period tt. As we mentioned above this is related to the nature of the TMS. Also this 
is manifested in the wavefunction ([5]) as linear Rabi frequency Q n in the parameter n. 
Moreover, from Figs. 2(a)-(b) the £(T) = at T = 0, and when the interaction is switched 
on, i.e. the interchange of energy between the two modes and atom occurs, the bipartite 
exhibits immediately maximum entanglement (i.e. £(T) provides maximum value), stays for 
a while then exhibiting long-lived disentanglement at T = tt/2. As the interaction proceeds, 
the £(T) goes gradually to the maximum value (i.e. entanglement) and eventually tends to 
zero at T = it. Furthermore, the comparison between Fig. 2(a) and Fig. 2(b) shows that the 
atomic coherence decreases the amount of entanglement in the bipartite. This agrees with 
the fact that the interference in phase space decreases the degree of entanglement in the 
bipartite j47[. To obtain more accurate information on the disentanglement in the bipartite 
for r) = we give a closer look at the analytical form of the wavefunction ([5]) at T = mir and 
T = m'7r/2, where m and m' are integer and odd integer numbers, respectively. At these 
values of the interaction times the wavefunction §5§ can be easily expressed as: 



\ip(T = mn)) = C n exp[irn(n+l)7r]|n, n)® 



n=0 



cos(-)| T) + sin(-)exp[ ? (0-m7r)]| j)j , (13) 
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\i){T = ^tt)) = £ exp(wwr) C^exp^) sin(f)|2n, 2n) 

n=0 L 



(14) 



-C 2n+1 exp(^)cos(f)|2n + l,2n + l) ® [exp(-z^)| T) + * -1 | I)], 

From ( jl~3l) it is evident that when T = mir the bipartite tends to its initial form, but the 
states of the field and atom include additional phase factors. This behavior is quite similar 
to that of the standard JCM [1] and two-atom JCM [47]. Nevertheless, from flUJ) the 
bipartite is disentangled, however, some information from the field is involved in the atomic 
system and vice versa. It seems that this is the origin in the long-lived disentanglement 
at these values of the interaction times. Moreover, expression ([14"]) shows that when the 
atom is initially in the excited (ground) state the system can generate odd (even) two-mode 
squeezed vacuum state. Thus at one-half revival time the system can generate two-mode 
cat states, which is similar to the standard JCM 0. We conclude this part by writing down 
the explicit form of the wavefunction for (8, 0, 77) = (vr/2, 0, 0), which can be evaluated as 

|V(T)) = -= C n exp[i(n + l)T)\n, n) ® [| T) + | 1)1 . (15) 

V 2 n=0 

From (1T51) it is evident that the bipartite exhibits always disentanglement, i.e. £(T) = 0. 
From the above investigation one can notice that for 77 = the behavior of the £(T) is 
consistent with that of the atomic inversion. 

Now we draw the attention to the influence of the Kerr-like medium on the behavior of 
the £(T) for the system under consideration (see Fig. 2(c)). The comparison between Fig. 
2(a) and Fig. 2(c) shows that involving Kerr-like medium in the system drastically changes 
the behavior of the £(T), e.g. it increases the amount of entanglement in the bipartite. 
This is related to several facts: (i) Kerr-like medium provides non-trivial photon-dependent 
phase in the wavefunction (cf. (j3J)) as well as it makes the Rabi oscillation as a nonlinear 
function in terms of n. (ii) The Kerr-like medium alone has a nonclassical nature where 
it can generate superposition states such as Yurke-Stoler states [48] . It is obvious that for 
77 7^ the evolutions of the (S Z (T)) and £(T) generally are not consistent. 

The final remark, for the system under consideration we have found that the von Neumann 
entropy gives typical behavior as that of the £(T). In the following section we study the 
evolution of the atomic Wehrl entropy and compare it with that of the £(T). 



IV. ATOMIC WEHRL DENSITY AND ENTROPY 

In this section we investigate the atomic Wehrl density and atomic Wehrl entropy AWE 
for obtaining more information on the system. Also we compare the behavior of these 
quantities with those given in section 3. We start the investigation by defining the atomic 
Q-function as 

Q a (Q, $, T) = ^ (0, $ \p a (T)\ 0, $) , (16) 

where p a (T) is the atomic reduced density matrix (jSj) and |0, $) is the atomic coherent state 
expressed as 

|0, $) = cos (0/2) |T) + sin (0/2) e i$ ||> , (17) 
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where < O < 7r, < $ < 2n. The definition (|T6|) means that two different spi n coherent 



states overlap unless they directed into two antipodal points on the sphere 1491. Here we 
shows that Q a can give qualitative information on the entanglement in the bipartite. From 
©, CEBD and (JT7J) one can evaluate Q a (0, $, T) as 



Q a (9,$,T) = — {l+p z (T) cos9 + [ Px (T) cos $ + p y (T)sm$) sine}. (18) 

One can easily check that the Q a given by ([TBI is normalized. From ([TBI it is obvious that 
at particular values of the G and $ the Q a can give information on the dipole moment and 
the atomic energy. Moreover, we have found that the behavior of the Q a {®, <&,T) confirms 
the information obtained in the previous section. For instance, when the atom is initially 
in the excited or ground state and T = mir, where m is integer, the Q a (Q,<&,T) exhibits 
behavior as that of the initial case. More illustratively, at these values of the interaction 
times we have p z (T) = p z (0), p x (T) = p y (T) = and hence (ITBl) reduces to 



Qa(6,$,T 



VTLTX] 



1 

47T 



[1 + p 2 (0)cos 0] 



(19) 



From (fT9~|) . Q a is independent of <3> and hence it exhibits wave shape in the G — $ plane 
with amplitude p 2 (0). On the other hand, when the atom is in the excited or ground or 
superposition states with (6,4>) = (tt/2,7t/2) and T = mix/2, where m is odd integer, the 
Qa(G, $, T) exhibits two-peak structure (, i.e. one peak is up and the other is down) as shown 
in Fig. 3a (for given values of the interaction parameters). Additionally, we plot Q a (Q, $, T) 
for the atomic superposition state at T = 7r in Fig. 3(b). The comparison between Fig. 
3(a) and Fig. 3(b) shows that the locations of the maximum and minimum values in the 
Q a (Q,$,T) are interchanged even though in these cases the bipartite is disentangled (cf. 
f fT4"|) ). For atomic trapped case (cf. (1151) ) we have p x (T) = p x (0), p y (T) = p z (T) = and 
hence (1181) becomes time independent. Now we draw the attention to the influence of the r\ 
on the behavior of the Q a (9, $,T), which is plotted in Fig. 3(c). We have generally found 
when the atom is in the excited or ground state and T = mir the Q a (Q, $, T) is insensitive of 
r), however, at T = rrnr/2, e.g. for the excited atomic state, the two-peak structure involved 
in the Q a for r] = is destroyed (compare Fig. 3(b) and Fig. 3(c)). This confirms that 
the bipartite is entangled (disentangled) at T = tc/2 (T = it). From the facts presented 
above and previous section we can conclude that Q a can give qualitative information on the 
entanglement in the bipartite. It is worth mentioning that the Q function of the field can 
manifest squeezing involved in the field components, however, this is not the case for Q a . 
This is related to the difference between the bosonic and fermonic algebras. 
Now we draw the attention to the AWE, which is defined as 33]: 



W a (T) = - 



2- 



Q a (0, $, T) In Q a (9, $, T) sin OdQd^. 



(20) 



It is worth reminding that the definition (1201) is given in analogy with that of the field Wehrl 
entropy (26[. Furthermore, the W a (T) cannot be negative as a result of the Q a is a non- 
negative function. As it is generally difficult to find a closed form for the W a (T) numerical 
techniques have to be used. Nevertheless, at particular values of the interaction parameters 
the exact form can be obtained. For instance, for f|T9l) the relation (1201) can be evaluated as 



W a (T 



mir 



p 2 M 



In 



i-p*(Q) 

_1 + P*(0)_ 



-In (1-^(0)). (21) 
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FIG. 3: The Q a (0,$,T) for sinh 2 r = 10 and (9,<f>,T,rj) = (0,0,vr/2,0) (a), (vr/2, vr/2, vr, 0) (b) 
and (0,0, vr/2, 0.03) (c). 
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FIG. 4: The W a (T) versus scaled time T for sinh 2 r = 10 and (9, cf), rj) = (0, 0, 0) (a) and (0, 0, 0.03) 
(b). 



From expression (12TT) for atom initially in the excited or ground state we obtain W a (T = 
m7r) ~ 2.488. In the language of entanglement the bipartite is disentangled when W a ~ 
2.488. This value represents the lower bound of the W a (T). In Figs. 4(a) and (b) we 
plot W a {T) for the same values of the interaction parameters as in Figs. 2(a) and (c), 
respectively. Surprisingly one can observe that the behaviors in the Figs. 4(a) and (b) are 
completely similar to those in the Figs. 2(a) and (c), respectively. We have checked this fact 
for different values of the interaction parameters and obtained the same conclusion. The 
similarity between the behaviors of the £(T) and W a (T) can be realized from the relations 
(fl2"|) and ({TBI) , i.e. the LE and the atomic Wehrl density and/or AWE are function in 
p x (T) , p y (T) , p z (T) . Moreover, the upper bound of the AWE can be obtained by comparing 
the maximum value of the LE (= 0.5) with the maximum value of the AWE (~ 2.53102). 
Thus the upper bound of the AWE is 2.53102. It is obvious that the AWE is compact 
informative one-parameter measure describing the entanglement dynamics of the system at 
all values of the interaction time. From this information we can conclude that the W a (T) is 
a new measure for quantifying the entanglement in the bipartite. 

The final remark is: the LE associated with the reduced atomic and reduced field density 
matrices can give typical information on the entanglement in the bipartite. This is not the 
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case for the field Wehrl entropy Wf and the atomic Wehrl entropy W a . For instance, as we 
mentioned in the Introduction the Wf gives information on the RCP occurred in the atomic 
inversion as well as on the splitting of the Husimi Q function in the course of the collapse 
region 3(j. The difference between the behaviors of the Wf and W a is related to that the 
Q functions of the field and atom provide different information on the system, where the 
former and the latter give information on the atomic inversion 50j and on the entanglement 
in the bipartite, respectively 



V. CONCLUSION 

In this paper we have studied the properties of the TJCM (involving a Kerr-like medium) 
when the modes and atom are initially in the TMS and atomic superposition state, respec- 
tively. Particularly, we have investigated the atomic inversion, LE, atomic Wehrl density 
and AWE. We have shown that all these quantities can give information on the entangle- 
ment in the bipartite. The bipartite exhibits instantaneous and long-lived disentanglement 
at particular values of the interaction parameters. We have derived the explicit forms for 
the atomic and field states at the disentanglement times. We have shown that the atomic 
Wehrl entropy (density) gives quantitative (qualitative) information on the entanglement 
in the bipartite. Most importantly, the LE and AWE can give typical information on the 
entanglement in the bipartite. This is related to that these quantities can be expressed 
in terms of the expectation values of the atomic operators. In this regard the W a can be 
interpreted as an information measure for such system. Finally, we have shown that the 
Kerr-like medium increases the degree of entanglement in the bipartite. 
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